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Abstract. Due to a new measurement of the S° mass, the Coleman-Glashow 
formula for the baryon octet e.m. masses (derived using unbroken flavor St/3) is 
satisfied to an extraordinary level of precision. The same unexpected precision 
exists for the Gell Mann-Okubo formula and for its octet-decuplet extension (G. 



Morpurgo, Phys. Rev. Lett. 68 (1992) 139). We show that the old question "why 



do they work so well?" is now answered by the general parametrization method. 



O 1 (PACS: 12.38.Aw; 13.40.Dk; 14.20.-c) 

a: 

1. Introduction 

A recent measurement of the S° mass |l] lowered considerably its error. The 

h ; 

H° mass is now 1314.82 ± 0.06 ± 0.2 MeV. The importance of a new measurement 
was noted long ago in connection with the Coleman-Glashow (CG) e.m. mass 
formula. The previous value was 1314.9 ±0.6 MeV [3]. Indeed now the agreement 
of CG with the data is more miracolous than ever. Writing the CG formula as: 

p-n = £ + -5r + (1) 

the present data give: 

l.h.s. = -1.29 MeV r.h.s. = -1.58 ± 0.25 MeV (2) 
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Because the mass differences X~ — S + in (|I]) is « 8 MeV, the agreement is amazing 
(= (4 ± 3)%) [before |l[], it was already excellent (1.29 to be compared with 
1.67 ± 0.6)]. To appreciate the point, note that the CG formula was derived 0] 
assuming unbroken flavor SU(3); but flavor is violated -in the baryon octet- by 
« 33%. 

A similar situation applies to the Gell Mann-Okubo mass formula and its octet- 
decuplet extension by one of us H. It also holds -with larger errors- for some 
formulas of Gal and Scheck ||. We already discussed all these relations using 
the QCD general parametrization method, but the result [p]] suggests a revisitation. 
Indeed we are dealing perhaps with one of the most precise estimates in processes 
where the strong interactions play a role. 

As stated above, the original derivation of CG neglected entirely the flavor 
breaking of the strong interactions. But it was shown in || that the CG formula 
can be derived also keeping all the flavor breaking terms, with the only omission 
of terms with 3-quark indices. Here we complete the derivation p[; we include, in 
addition to the terms considered in 0, the effect of the m d — m u mass difference 
and the so called Trace terms, absent in ; they do not alter the conclusions of |2| . 

2. A brief summary of the general parametrization method 

It is convenient to recall briefly the QCD parametrization method |7|, ||. The 
method, based only on general features of QCD, applies to a variety of QCD ma- 
trix elements or expectation values. By integrating on all internal qq and gluon 
lines, the method parametrizes exactly such matrix elements. Thus hadron prop- 
erties -like e.m. masses, including their flavor-breaking contributions- are written 
exactly as a sum of some spin- flavor structures each multiplied by a coefficient. 
Each structure (term) has, for baryons, a maximum of three indices. The coef- 
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ficients of the various terms are seen to decrease with increasing complexity of 
the term. By the way this "hierarchy" explains why the non relativistic quark 
model (NRQM), that keeps only the simplest terms, works quantitatively fairly 
well. Though the parametrization is performed in a given Lorentz frame and is, 
therefore non covariant, it is relativistic, being derived exactly from a relativistic 
field theory, QCD. For the basis of the method see |7|, applications are also 
given in [2,5,7-11]. Other references are listed in (JTTJ, |T2|; the latter gives a short 
review. 

Here we will not recall the details of the method, but -for completeness- sum- 
marize it. The e.m. contribution to the mass of a baryon B is: 

(Vb|J#b) = (<Pb\V^V\<Pb) (3) 

In (3) Q is -to second order in the charge- the exact QCD operator, expressed in 
terms of quark fields, representing the e.m. contribution to the mass, including all 
flavor breaking contributions of the strong interactions; \i[>b) is the exact eigenstate 
of B at rest of the QCD Hamiltonian; \4>b) is an auxiliary three body state of B, 
factorizable as 

\4>b) = \X L =o ■ W B ) (4) 

into a space part Xl =0 with orbital angular momentum zero and a spin unitary- 
spin part Wb- The unitary transformation V -applied to the auxiliary state \4>b)- 
transforms the latter into After integration on the space variables, (4) can 

be written 

{M^b) = (MV^V^b) = {Wb^UT^sJWb) (5) 

V 

where T v are operators depending only on the spin and flavor variables of the 
three quarks in 4>b and the t u 's are a set of parameters. Of course the operator 
V dresses the auxiliary state \4>b) with qq pairs and gluons and also introduces 
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configuration mixing. Thus it generates the exact QCD eigenstate \i[>b)- It is the 
factorizability of \<$b) that allows the second step in Eq.(5), eliminating the space 
coordinates. The list of r^(s, /)'s in (5) for the e.m. masses was given in 0. As 
stated, we will revisit this parametrization. 

3. The corrections to the octet masses in the Coleman-Glashow formula 

The e.m. corrections to the octet baryon masses including flavor breaking 
analyzed in left out the rrtd — m u contribution and the Trace terms that, in this 
paper, we will include later (Sect. 4). For the moment we reanalyze the results of 
||. There we called SiB the e.m. mass correction of baryon B (i — 0,1,2 refers 
to the order of the s-flavor breaking). Thus by 5qB we meant (and mean) the 
e.m. correction neglecting all flavor breaking effects, S±B is the e.m. correction 
including only first order flavor breaking effects and so on for 82B. 

The quantities T u (s, /) that enter in the construction of <5;.B's are given in Eqs. 
(17-19) of [Q. There is nothing to change in these equations except for a (trivial) 
point of notation. In |2| the strange quark was called A (the non strange ones M 
and V). Here we use the standard current notation s, d, u. We do this because 
(see H) we can select as we like the q 2 of the renormalization point of the quark 
mass and we now turn to the standard q around 1 GeV. Thus the projectors P x , 
P M , P v in @ are now written P s , P d , P u . 

The IVs in Eqs. (6,7) below are the same as those in Eqs. (17, 18) of (the 
sum symbol in each T u is defined in ||); we transcribe them here: 

1) T's of zero order in flavor breaking: 

ri = E[Q?] ; r 2 = E[Q 2 M • ; r 8 = £[<#fo ■ <r*)] (6) 



4 



2) r's of first order in P s (acting in A,S,E*,H,H*,0): 

r 7 = TXQlPi] ; r 8 = £[Q?jy(<r, • <r*)] ; r 9 = E[Q?/7fo • 

rio = E[QM ; r u = ElQ^to • <r fc )] ; r 12 = E[Q?nV* + **) • os)\ 

ri3 = E[Q<Q*^ 8 ] ; r 14 = E[QiQ*J?to • ; r 15 = EfQiQ^K + ^) • ^1 

r 16 = E[<5iQitP/] ; r 17 = EiQiQkPfoi ■ <r*)] ; r 18 = EtQ^/K + <r k ) ■ 

(7) 

where Qi = ~P" — |P/ — ~P/ and the Pj's are the projectors on the u, d, s quarks. 

As to the IVs of second and third order in P s , they are listed in Eqs. (19,20) 
of 0. We will not transcribe them here. 

The expression 5qB of the electromagnetic mass of B at zero order in flavor 
breaking is: 

5 q b = oTi + 6r 2 + c r 3 + dr 4 + e r 5 + /r 6 (8) 

(where, to agree with p|, we used a . . . f instead of t\ . . . t§). As shown in p| one 
can check that: 

5 P - 6 N = 5 S+ - 5 £~ + 6 ~~ - 6 ~° 

which is the CG relation at zero order in flavor breaking. 

As to 8\B and S2B, coming from the first and second order in flavor breaking, 
we refer to 0. There it is shown (table III) that, except for terms with three 
indices, also 5\B and <5 2 P leave the CG relation unaltered. Thus, to all orders in 
flavor breaking, with the only omission of three quark terms in 5\B and <5 2 P, the 
CG relation holds: 

5P — 5N = <5£ + - 6E- + 8E~ - 5E° (9) 

where now 5 = 5o + Si + 5 2 . 

To evaluate the order of magnitude of the three quark terms of the type 5\B 
and 82B (the only ones that violate the CG formula) we now use the hierarchy 
discussed at length in previous papers ||, |8|. Of course, the dominant order of 
magnitude, is here that of the two index terms of the form QiQk- Let us consider 
the three quark terms in 8\B. One can show that r 9 , r i2 and Ti 5 do not contribute 
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to the left and right hand sides of the CG formula. As to r 16 , r 17 and r 18 , they do 
not contribute to 8\n and S\p, while for — <5iX _ + 5iE~ — <5iH° one gets in all 
cases ("something" )/3. As to the magnitude of "something" , note that according 
to the hierarchy, each of Ti6, Tu and Tig carries a reduction factor « 0.3 due to 
the presence of a P s and w 0.3 due to one more gluon exchange. Thus we have 
a reduction of order (1/3) 3 = 4 ■ 10~ 2 for each of the above terms with respect to 
the dominant no flavor-breaking contributions. (Note: r 18 here is 1/2 that in 0] 
due to an incorrect normalization in p that produced, however, no effect because 
3-index terms were not evaluated in ||.) 

Because experimentally S~ — S + = 8 MeV and 2~ — 5° = 6.4 MeV, the above 
value of 4 ■ 10 -2 implies an expected difference between left and right hand sides 
of the CG formula of ~ 0.2 0.3 MeV that does not disagree with the data. 

A similar argument holds for the three-quark terms of second order in flavor 
breaking listed as (19) in 0. These second order flavor breaking three-quark terms, 
are expected from the hierarchy, to be ~ (1/9) of the first order flavor breaking 
terms mentioned above; that is 

\5 2 E- - 5 2 ~°\ « l\5^ + - <JiE" + 5iS- - <fiS°| 

and therefore the three quark, second order flavor breaking effect should contribute 
to the difference between left and right hand side of the CG relation by ~ 0.02 
MeV. 

4. The effects of m^ — m u and of the Trace terms on the CG relation 

a) The effect of m^ — m u . As shown in || the quantity that intervenes in 
evaluating the effect of a mass difference Am between quarks is Am/(j3A) where 
A = Aqcd — 200 MeV and (3 is approximately 3. Because Am for d and u is a 
few MeV, only the first order term in the expansion in Am/ (13 A) may be relevant. 
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But it is easy to check that the CG formula is not modified. (Obviously product 
terms of Am and QiQk type perturbations are totally negligible.) 

b) The effect of the Trace terms. In addition to the terms in M other terms are 
present in the general parametrization (see in particular footnote 14). They 
leave the CG formula unaltered, as we will show; they must however be recorded. 
These "Trace" terms correspond to QCD Feynman closed loops, as exemplified, 
e.g. in [TIJ (fig.l) or |T2| (fig.3). Introducing the matrix 

Q = -P u - -P d - -P s (10) 
3 3 3 

(not to be confused with the baryon charge Qb that was called Q in [g]), the Trace 
terms can be constructed as follows: Consider the quantities 

T x = Tr{QP s ) , T 2 = Tr(Q 2 ) , T 3 = Tr{Q 2 P s ) 

and multiply them by J^i&i ■ CT fc) or by the sums listed in Eqs.(6) and (7) keeping 
only those expressions that, after the multiplication, contain two Q symbols (either 
Q 2 or Q ■ Qi or QiQk) and a number of P s from to 2; for instance (just to 
exemplify) 

Tr(Q 2 ) E(cTi ■ cr k ) , Tr(Q 2 P s ) . <j k ) , Tr(QP s ) EQi , n n 
Tr(QP s )j:Q t (^-<Tk) , Tr(QP s )j:QA^-^k) 1 ' 

It is easy to check that none of the above terms changes the previous conclusions 
concerning the exactness of the CG equation. The only new quantity that en- 
ters produced by terms of type (11) is the expectation value of J2i Qi&zi on the 
baryons. For instance Tr{QP s ) J^Qii&i ' °"fe) — —Qb/2 + (3/2) Y^iQi^zi] the val- 
ues of (J2iQi°~zi) are listed in Table I of 0, indicated there as (E|). One checks 
immediately that the CG equation remains true. 

We conclude this discussion of the CG miracle as follows: The three quark 
terms are expected to give a very small, but non zero contribution to CG. It is 
remarkable that the hierarchy typical of the general parametrization, appears to 
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explain this smallness thus providing one of those few cases where one can estimate 
a tiny effect of the strong interactions and find it compatible with the data. 

5. The Gell-Mann Okubo mass formula and its extension to the octet- 
decuplet 

We will discuss now, first qualitatively, then quantitatively, the reason why, 
besides the CG formula, also the Gell-Mann Okubo (GMO) formula has its share 
of mysterious perfection. The GMO formula for the octet baryons: 

i(n + S°) = l(3A + S) (12) 

is derived neglecting terms of second order in flavor breaking. The expansion 
parameter for flavor breaking is -as generally known and determined from the 
general parametrization 0, ||- 0.3 to 0.33. Because the magnitude of first order 
flavor breaking is ~ 150 to 190 MeV, an estimate of the contribution of second 
order flavor breaking terms -neglected in (12)- is ~ 50 to 60 MeV. Instead the r.h.s 
and l.h.s of (12) differ by ~ 6 MeV. Is this success just luck? Not entirely. In fact, 
by the general parametrization, we can now parametrize together the decuplet and 
octet masses and determine in this way all the 8 parameters in the octet+decuplet 
mass formula. The coefficient of the second order flavor breaking, (a + b) in f|], is 
in fact, as expected, ~ 3 times smaller then the above estimate 50 or 60 MeV, due 
to the fact that necessarily second order flavor breaking terms have two indices. 
Thus one expects a difference of = 17 or 20 MeV between the left and right hand 
sides of (12). Because || the parameter that enters in the GMO formula is (a+b)/2 
we are led, by this order of magnitude argument, to a difference of 8 — 10 MeV, 
not far from the experimental value of ~ 6 MeV. The essential role is once more 
played by the hierarchy, that again leads to negligible three index second order 
flavor breaking terms, precisely as for the CG formula. 
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Let us now be more quantitative. The fact that the coefficient c, d (Eq. (4,5) 
of ||) multiplying the second order flavor three index terms are negligibly small 
was first discovered in relation to the GMO formula (see also [Q) and used 
above in discussing the CG formula. Barring these c, d coefficients one finds a 
mass formula that relates the octet and decuplet masses, correct to second order 
in flavor, except for a three index term. This formula is just || the GMO formula 
plus a "decuplet" correction T: 

hn + S°) + T = i(3A + 2E+ - S°) (13) 

where S + = 2S + — S° + 2(p — n) and T is: 

T = H— - -(ft + E*~) (14) 

The charge specification are inserted here because, at this accuracy, Eq.(13) must 
take into account the e.m. contributions: The combination (13) is unaffected by 
the e.m. corrections if the latter are calculated at zero order flavor breaking. (Note: 
Eq.(13) is more simply |(p + S°) + T = |(3A + 2S+ - S°) but here we kept the 
form used in 0.) 

While of course the improved S° []J slightly decreases the experimental error 
in (13), the main part of the error comes from the decuplet masses in T. With the 
conventional values of the masses [§] it is: 

T [conventional) = 5.18 ± 0.66 MeV 

whereas, if the pole values || of the resonances are taken ||, it is: 

T(pole) = 6.67 ± 1.25 MeV 

The left and right hand side of (13) become: 

(conventional) Left = 1132.36 ± 0.7 MeV Right = 1133.93 ± 0.04 MeV 
(pole) Left = 1133.86 ± 1.25 MeV Right = 1133.93 ± 0.04 MeV 
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In both cases it is again true that three quark terms breaking flavor to second 
order are estimated to contribute to the difference between r.h.s. and l.h.s. less 
than 0.7 MeV. 

We finally note that, to first order in \m u — m^j (/3Aqcd), the u, d mass 
difference does not affect the octet-decuplet mass formula (13); also the Trace 
terms do not modify the parametrization in this case. 

6. Other e.m. mass formulas 

Using the NRQM Gal and Scheck |J derived long ago a set of relations among 
the electromagnetic masses for mesons and baryons. For mesons their assumptions 
were very restrictive, but for baryons they amounted mostly to the neglect of three 
body terms. In || we showed how these formulas could be reproduced by the 
general parametrization method under certain assumptions. Although it might be 
of interest to reanalyze the situation in more detail, we refrain from it because the 
experimental data did not change substantially. 

However, in order to stimulate more precise measurements (if possible), we 
write down below a relation that is totally analogous, for the decuplet, to the 
Coleman Glashow formula for the octet and that can be treated in exactly the 
same way. It is independent of the Gal Scheck derivation, but easily verifiable 
using the Eqs.(27, 28) of [§: 

5A + - 5A° = - 8E*~ + 8E*~ - 5E*° (15) 

Incidentally Eq.(15), together with the relation (16) below, derived, as well known, 
from isospin algebra 

SA++ - 5 A- = 3(5 A + - 8 A ) (16) 
might be of interest in connection with the determination of the mass differences 
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between the A's. More generally, Eq.(15) can be derived, using (16), from the 
charge corrected second order Okubo equations (7, 8) of ref . IfUfl . 



7. Conclusion 

The reason for the extraordinary perfection of the Coleman- Glashow relation 
that holds at present to = (0.29 ± 0.25) /8 = (4 ± 3) • 10~ 2 (in spite of having been 
derived, we recall, using unbroken flavor SU(3)) is now clear. It depends, we have 
shown, on the smallness of the three index terms in the general parametrization. 
We underline again (Sect. 4) that this smallness represents one of the few cases 
where, thanks to the hierarchy in the parametrization, an estimate of an effect due 
to the strong interaction can be given and found to be tiny as expected. 

As to the Gell-Mann Okubo formula, its octet decuplet extension including 
second order flavor breaking except for three index terms, holds to better than 
2 • 10~ 2 . This value is the ratio between the experimental error (« 1 MeV) and 
what one would expect estimating the orders of magnitude (~ 50 MeV) just by 
flavor breaking to second order. This confirms the smallness of the three index 
terms. 
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